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1 Introduction 



Integrable quantum systems containing Fermi fields have been attracting increasing interest due to their 
potential applications in condensed matter physics. The prototypical examples of such systems are the 
supersymmetric generalizations of the Hubbard and t-J models [1]. They lead to a generalization of the 
Yang-Baxter (YB) equation [2] associated with the introduction of the a Z2 grading [3]. In addition to 
the grading, it is also important to introduce open boundary conditions to study the boundary effects 
on the bulk system. The most powerful method in the analysis of integrable models is the Bethe Ansatz 
(BA). The algebraic BA, also known as the quantum inverse scattering method (QISM) [4], is an elegant 
and important generalization of the coordinate BA [5] 

In this work wc will study two graded three-state vertex models with reflecting boundary conditions. 
Their boundary algebraic BA are delineated based in the recent progresses [6, 7], for the formulation of 
the QISM for the 19-vertex models with boundary condition terms. 

Let F = Vo© Vi be a Z2-graded vector space where and 1 denote the even and odd parts respectively. 

s 

Multiplication rules in the graded tensor product space V ^ V differ from the ordinary ones by the 

s s 

appearance of additional signs. The components of a linear operator B G V result in matrix 
elements of the form 

{A 4 B fJ^ = (-)K/3)(p(a)+p(7)) A^^Bps. (1.1) 

S S 

The action of the graded permutation operator V on the vector |a) (8) eV<S)Vis defined by 

P \a) ® 1/3) = {-)pMpW 1/3) ^ \a) =^ {VfJ^ = {-f^^^^^^H^s 5^^. (1.2) 
The graded transposition st and the graded trace str are defined by 

(^'*)a/3 = (-)(^(") + ^)^(«^/3a, Str^ = (1.3) 

a 

where p{a) — 1 (0) if \a) is an odd (even) element. 
For the graded case the YB equation 

ni2{u)nii{u + v)n2s{v) = 7^23(v)7^l3(« + v)ni2{u) (i.4) 

and the reflection equation [8, 9] 

7^l2(^^ - v)K^{u)TZ2i{u + v)K2{v) = Jfa" (^;)7^l2 (« + v)K^{u)n2i{u - v) (1.5) 

remain the same as in the non-graded cases and we only need to change the usual tensor product to the 
graded tensor product. 
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In general, the dual reflection equation which depends on the unitarity and cross-unitarity relations 
of the 7?,-niatrix takes different forms for different models. For the models considered in this paper, we 
write the graded dual reflection equation in the form introduced by Zhou at a/. [10] (see also [11]): 

ni\''''{-u + v){K+y*'{u)M^^nil''''{-u -v- 2p)Mi{K+y*'{v) 

= {K+y^^{v)Mxnf^'*^{-u -V- 2p)M^\K+y*^ (u)7e^*i^*^(-M + v), (1.6) 

and we will choose a common parity assignment: p(l) = p(3) — and p{2) = 1 , the BFB grading. 
Now, using the relations 

nil^'^^u) = hR2i{u)Ii, 7efi^'*^(«) = /ii?i2(u)/i and IK+{u)I = K+{u) (1.7) 

with I = diag(l, —1, 1) and the property [MiM2,Ti{u)] = we can see that the usual isomorphism [12] 

K-{u) :^ K+{u) = K-{-u - py*M. (1.8) 

holds with the BFB grading. Here stj denotes super-transposition in the space i. 

A quantum-integrable system is characterized by the monodromy matrix T{u) satisfying the funda- 
mental relation 

R{u - v) [T{u) ^ T{v)] = [T{v) (8) T{u)] R{u - v) (1.9) 

where R{u) is given by R{u) = PTZ{u). 

In the framework of the QISM [4], the simplest monodromies have become known as C operators, the 
Lax operators, here defined by Caq{u) = TZaq{u), where the subscript a represents the auxiliary space, 
and q represents the quantum space. The monodromy matrix T{u) is defined as the matrix product of 
N Lax operators on all sites of the lattice, 

T{u) = CaN{u)CaN-l{u) ■ ■ ■ Ca\{u). (1-10) 

The main result for open boundaries integrability is: if the boundary equations are satisfied, then the 
Sklyanin's transfer matrix [8] 

t{u) = str^ {K+{u)T{u)K-{u)T-^{-u)) , (1.11) 
forms a commuting collection of operators in the quantum space 

[t{u),t{v)]=(), yu,v (1.12) 



2 



The comniutativity of t(u) can be proved by using the unitarity and crossing-unitarity relations, the 
reflection equation and the dual reflection equation. In particular, it implies the integrability of an open 
quantum spin chain whose Hamiltonian (with K~{0) = 1) is given by [8] 

H- - V" H- -L^ dK-{u) 

fc=l "=o 
where the two-site terms are given by 



stroj^+(0)gjv,o 
strii:+(0) ' 



(1.13) 



(1.14) 



u=0 



in the standard fashion. 

The paper is organized as follows: in Section 2 we define the models to be studied. In section 3, 
we present our detailed calculations common to both models and in Section 4 the eigenspectra and 
the corresponding Bethe equations are explicitly presented for each model. Section 5 is reserved for 
conclusions. 

2 The models 

The three-state vertex models that we will consider are the sZ(2|l)(^) model and the osp{l\2) model. 
Their 7^-matrices have a common form 

\ 





/ Xi 


X2 


X3 


X5 


xe 




X7 


■R{u) = 




2/5 


ye 


X2 




X2 


X5 




\ 




2/7 




ye 


2/5 


X3 

X2 

XI 



(2.1) 



satisfying the properties 



regularity 
unitarity 
PT — symmetry 
crossing — symmetry 



7ei2(o) = /(o)i/2p,2, 
■R,2{u)nir'H-u) = f{u), 

Pl27^12(w)Pl2=7^f2^^*M«), 

n^2{u) = UiTZti{-u-p)U^\ 



(2.2) 



where f{u) = Xi(u)xi(—u), Xi(u) being defined for each model below, p is the crossing parameter and 
U determines the crossing matrix 



M = WU = M* 



(2.3) 
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Unitarity and crossing-symmetry together imply the useful relation 

Mi<2^ {-u - p)M^^nll' {u-p) = f {u). (2.4) 
2.1 The sl(2|l)(2) model 

The solution of the graded YB equation corresponding to s/(2|l)(^) in the fundamental representation 
has the form (2.1) with non-zero entries [13, 14]: 

x\{u) = sinh(tt + 2?7) cosh(M + 77), a;2(M) = sinhu cosh(it + 77), 

X3{u) = sinhucosh(tt — 77), a:4(tt) = sinhucosh(tt + 77) — sinh2r7COshry. 

2/5(u) = a:;5(M) = sinh 277Cosh(u + 77), 2/6(u) = a;6(M) = sinh27/sinhu, 

y7{u) = a;7(M) = sinh2r/cosh77. (2.5) 

This 7^-matrix is regular and unitary, with f{u) = xi{u)xi{—u), P- and T-symmetric and crossing- 
symmetric with M = 1 and p = f]. The graded version of the crossing-unitarity relation (2.4) is satisfied 
with /'(u) = xi{u + i^)xi{—u — i^). 

The most general diagonal solution for K~{u) has been presented in Ref. [15] and it is given by 



K-{u,P)={ 1 I, (2.6) 



^33(") 



with 



/3 sinh u + 2 cosh u /3 cosh(u + 77) — 2 sinh(u + 77) 

ii^"'' = ~/3sinhw-2cosh«' ''^al^) = ^ ^^^^^ _ ^) ^ 3 sinh(w - 7?) ' ^ ' 



where /3 is a free parameter. Due to the automorphism (1.8) the solution for K^{u) is given by K {—u- 
p, \a) I.e. 

K+{u,a)=\ 1 1, (2.5 



where 



r \ acosh(M + 7;) — 2 sinh(M + ry) ]^+ ( ) a sinh m + 2 cosh u 

(."J - cogii^y + r]) + 2 sinh(u + 77) ' 33W - q, sinh(u + 27?) - 2 cosh(u + 2??) ' ^39) 

and a is another free parameter. 
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2.2 The osp(2|l) model 

The trigonometric solution of the graded YB equation corresponding to osp{l\2) in the fundamental 
representation has the form (2.1) with non-zero entries [14]: 

Xi{u) = sinh(u + 2r]) sinh(u + Sr]) , X2 (w) = sinh u sinh(u + 3?7) 

xs {u) = sinh u sinh('u + ??) , X4 (u) = sinh u sinh(M + 877) — sinh 2r] sinh 3r] 

X5{u) = sinh 2?7 sinh(M + 877), 2/5 (m) = e" sinh 277 sinh(u + Sry) 

xq{u) = —e""^^'' sinh 2?7 sinh ?i, (u) = e"'*'^'' sinh 2?7 sinh w 

X7{u) = e"'' sinh 277 (sinh('u + 877) + e^'' sinh u) 

t/7(u) = e" sinh 277 (sinh(u + 877) +e'' sinh u) (2-10) 

This 7?.-matrix is regular and unitary, with f'{u) = f{u) = Xi{u)xi{—u). It is PT-symmetric and 
crossing-symmetric, with p = 877 and 



e 



-2r, 



M = \ 1 . (2.11) 

e^" 

Diagonal solutions for K~{u) have been obtained in [16]. It turns out that there are three solutions 
without free parameters, being K~{u) = 1, K~{u) = and K~{u) = F~, with 

/ ^e-2«/W(^) \ 
F^=\ 1 , (2.12) 

where we have defined 

= + ;(-)(.) = cosh(7. + 3,/2) 

smh(?i — 8r//2) cosh(u — 8?7/2) 

By the automorphism (1.8), three solutions K~^{u) follow as K'^{u) = M, K'^{u) = G+ and K~^{u) = G~ , 
with 

/ ^e2«+4'75(±)(w) \ 
G± = 1 , (2.14) 

where we have defined 

^ ^"^-sinh(u + 977/2)' ^ ^"^-cosh(u + 977/2)- ^^"^^^ 

We have thus nine possibilities for the commuting transfer matrix (1.11). We will only consider 
three types of boundary solutions, one for each pair {K~{u), K~^{u)) defined by the automorphism (1.8): 
(1,M), (F+,G+)and {F-,G-). 
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3 Algebraic Bethe Ansatz 

The monodromy matrix T{u) (1.10) and its reflection T~^{—u) can be written as matrices 3 by 3 

/ Tn(u) TM TM \ ( Tr,\-u) T^^\-u) T,f{-u) 

T{u)=\ TM T^siu) , T-\-u)=\ T^\\-u) T,-,\-u) T^-,'{-u) 

V niiu) Ts2{u) T^siu) J \ T,-,\-u) T,^\-u) T,,\-u) J (3.1) 

where 

T,,{u) = Yl 4?K^) ^ 4f.;'^K^) ^ • • • ® 4L,aM (3.2) 

fel,...,feiV-l = l 

where C\^^ are 3x3 matrices acting on the nth site of the lattice, defined by 

xi 

£(7) = I X2 I , C^^^={ X, I , 
0:3 




?/5 

£(7) = I J/6 I , C^^={ X4 1 , 




-31 ~ 








2/7 \ 


' -^32 ~ 1 


^ 


2/6 


° 1 


/-(») _ ( 
' ■^33 ~ 1 

















\ 





/ 




V 





/ 






(3.3) 



Using the unitary relation in (2.2) we can see that the reflected monodromy T ^{—u) has the following 
matrix elements 

Ttfi-^) = JT^ E <l\ i-u, -V) ® 41(-"> -r?) ® • • • 4 cZl^h^' -^)- 

•^W fc,,...,fc„_,=i (3.4) 

For the vertex models considered in this paper we can choose the highest weight vector of the mon- 
odromy matrix in a lattice of A'' sites as the even (bosonic) completely unoccupied state 

|0) = n®|0)fe, |0),= (oy (3.5) 
fc=i \ / 

where |0)j, is the local reference state at the k-th lattice site with three components 
The action of T{u) and T~^{—u) on this state are 

Xi{u) |0) * ** 

T(«)|0) = /^(«)T-i(-«)|0)= ( x^{u)\0) *** |, (3.6) 

x^{u) |0) 
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which give us, in the usual BA language, creation and annihilation operators for (3.5). Moreover, we are 

working out with boundaries and in this case we have a double-row monodromy defined by 

/ Uii{u) Ui2{u) Uu{u) \ 
U{u)=T{u)K-{u)T-\-u)= I U2i{u) U22{u) C/23H (3-7) 

V Usiiu) U32{U) U33{u) J 

where K'^^^u) is a reflection matrix. 

For K^~\u) = diag(A:fi(w), fc^2(")i ^33(^))' the matrix elements of U have the form 

3 

Uij{u) = J2Tia{u)kMT-j'{-u), i,i = 1,2,3. (3.8) 

a=l 

It follows from (3.8) that we will need to know the commutation relations between the operators T{u) 
and T~^{—u) in order to get the action of U{u) on the reference state (3.5). Using the fundamental 
relation (1.9) with u = —v we will get the matrix relation 

T^\-u)Ri2{2u)Ti{u) = Ti{u)Ri2{2u)T2\-u) (3.9) 



Applying both sides of this relation on the reference state, we find the following relations 



T2,{u)T^2\-n) |0) = /,(u)^^W^^ |0) (3.10) 



Ts,{u)Trs\-u) |0) = |/2(«)^^ - M^)hi^)jW^ - /^W^v^l |0) (3.11) 

u) — X', 

1^ 



n2{u)T,-,\-u) |0) = f2iu)^^j^j^^ |0) (3.12) 



where 



f I N 2/5 (2w) , , . Vri^u) 

hiu) = —TTT-:, f2[u) - 



fsiu) = 



xi{2u)y^{2u) — xr,(2u)y-j{2u) 
xi{2u)x4{2u) + X5(2u)y5{2u) ' 



^ X4{2u)y7i2u) + yU2u) 
' xi{2u)xi{2u) + xn{2u)yr,{2uy ^' ' 

Using these relations we can get the action of each operator Uij on the reference state: for the diagonal 
entries we have 

^ xl^ju) 



Un{u)\0) = A;n(«)^^|0) 



U22{u) |0) = /i(u)[/n(«) 10) + [k^^{u) - KMhiu)] |0) 



U33{u)\0) = [{f2{u)-.fi{u)Mu))Un{u) + f3{u)U22{u)]\0) 

+ few - kMMu) - fcn(«)/4(«)] |0) (3.14) 
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and for the elements out of the diagonal we get 

Uij{u)\0)=0, {i>j), C/i,(u)|0)^{0,|0)}, {i<j) (3.15) 
Now, we define news operators: 

V^{u) = U^^{u), Bi{u) = Ui2{u), B2{u) = Uis{u) 

=t/2l(u), V2{U)=U22{U)- fl{u)Vi{u), S3 («) = t/23 («) 

C2{u)=U3liu), C3{U)=U32{U), I?3 («) = C/33 (w) " /2 («) - /a (u)P2 (u) 

(3.16) 

to write the double-row monodromy matrix in the form 

/ Vi{u) Bi{u) B2{u) 
U{u)^U{u)=\ Ci{u) V2{u) B^iu) |. (3.17) 
V C2{u) C3(w) V3{u) 

The action of U (u) on the reference state has the usual BA form 

A'i(u)|0) * ** 
U{u)\Q)=\ X2{u)\Q) *** I (3.18) 
X3{u) |0) 



where 



X2{u) = [A;22(w) - A;ii(w)/i(u)] 



W = [k,,{u)~k^,{u)Mu)-k^,{u)U{u)]jj^ (3.19) 

The transfer matrix t{u) = str{K~^U), with diagonal left reflection K^~^^ = diag(A;^j, Aij^.fcg^) and 
BFB grading, has the form 

t{u) = kti{u)Un{u) - k+2{u)U22{u) + k+^{u)U33{u) 

= ni{u)Vi{u) + n2{u)V2{u) + n3{u)V3{u) (3.20) 



where 



^liu) = k^iiu) - /i(M)fc+ (m) + /2(u)A;^3(-«) 
^2{u) = -kt2{u) + f3{u)k^s{u) 

Qsiu) = fc+3(u) (3.21) 
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Here we note the sign (— ) from the graded trace is absorbed in the definition of Cl2{u)- 

From U{u) foUows the usual algebraic BA structure. Therefore we can look for states created by the 
operators Bi{u) from a reference ^"0 which will be eigenstates of (3.20). To do this we first recall the 
magnon number operator 

JV 

M = 5^ Mfe, Mk = diag(0, 1, 2) (3.22) 

k=l 

This is the analogue of the operator used in the coordinate BA construction. The relation m = N — S^ 
(M^'m = m'^ni ) allows us build states ^'m such that t{u)'^rn = A^^'m. Therefore, we can start the 
diagonalization of t{u) by considering all possible values of m in a lattice with sites. 

By the previous construction, "ifo is our reference state |0 > , which is itself an eigenstate of t{u) 

t{u)^o = Ao(w)*o (3.23) 

with eigenvalue 



Ao(u) = [ktiiu) - h{u)kt2{u) + /2(w)A;^3(w)] k-^-^{u) 



+ [-^iW + h{u)kts{u)\ [k22{u) - A;ii(m)/i(m)] 



+k+iu) few - kMMu) - k^,{u)U{u)] (3.24) 
It is the only state with m = 0. 

3.1 The one-particle state 

For m = 1 we seek a state of the form 

*i(ui)=Bi(ui)|0). (3.25) 
The action of t{u) on this state is given by 

t{u)^i{ui) = fli{u)Vi{u)Bi{ui) |0) + fl2{u)V2{u)Bi{ui) |0) + 03(^)2^3 (w)Bi(ui) |0) . 

(3.26) 

Since we know from (3.18) the action of the operators X>i(u) on the reference state |0), we need to arrange 
the operators products 

Pi(w)Si(wi), I?2(w)Si(wi) and I?3(w)Bi(wi) (3.27) 

in a normal-ordered form [17]: We anticipate that, in general, the operator- valued function ^'„(ui, . . . ,u„) 
for a n-particle Bethe state will be composed by a set of normal-ordered monomials. A monomial is said 
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to be in normal order if all elements Bi are on the left, and all elements Ci are on the right of the elements 
■Di. 

In order to get this normal ordering we recall that the double-row monodromy matrix satisfies 
the fundamental reflection equation 

■JZi2{u-v)Ui{u)n2i{u + v)U2{v) =U2{v)nx2{u + v)Ui(u)n2i{u-v), (3.28) 

where Ux{u) = U{u) (g) 1, U2{u) = 10 U{u) and '7^2l(u) = VTZi2{u)V. 

In the appendix of [7] it was shown as this equation can be used to recast the non-normal ordered 
products as the above into a linear combination of normal-ordered ones. However, in this article we will 
present only the main results. (See below) 

For the present case f(u)^'i(ui) can be computed with the aid of the following normal-ordered (or 
commutation) relations 

I>i(m)Si(mi) = aii{u,ui)Bi{ui)Vi{u) + ai2(u, mi)Si(u)I?i(ui) + ai3(u, ui)Si(u)I?2(ui) 
-|-ai4(u, ui)B2{u)Ci{ui) + ai5{u, ui)B2{u)C3{ui) + aie{u, ui)B2{ui)Ci{u) 

(3.29) 

T>2{u)Bi{ui) = a2i{u,ui)Biiui)T>2{u) + a22iu,ui)Biiu)Viiui) + a2a{u,ui)Biiu)D2iui) 
+a24{u,ui)B3{u)T)i{ui) + a25{u,ui)B3{u)T)2iui) + a2e{u,ui)B2{u)Ci{ui) 

+a27{u, Ui)B2 (u)C3(ui) + 028(W, Ui)B2 (Wl)Ci (u) + 029 («, Ui)B2{ui)Cs (u) 

(3.30) 

I?3(u)Si(ui) = a3i(u, ui)Bi{ui)'Ds{u) + as2{u, ui)Si(u)X>i(ui) + as3{u, ui)Si(u)X>2(wi) 
+a34{u,ui)B3{u)'Di{ui) + a35{u,ui)B3{u)'D2{ui) + a3e{u,ui)B2{u)Ci{ui) 
+a3r{u,ui)B2{u)C3{ui) + a3s{u,ui)B2{ui)Ci{u) + a3g{u,ui)B2{ui)C3{u) 

(3.31) 

Substituting these relations into (3.26) one gets 

t{u)^i{ui) = fli{u)Vi{u)Bi{ui) \0) + fl2iu)V2{u)Bi{ui) \0) + fl3{u)V3{u)Bi{ui) |0) 

= [an (u, ui)Oi {u)Xi {u) + a2i [u, ui)Q.2 {u)X2 [u) + 031 (u, ui)Vl3 {u)X3 {u)] *i (ui) 

3 3 
+ [Xi{ui)'^Vtj{u)aj2{u,ui) + X2{ui)'Y^D.j{u)aj3{u,ui)]Bi{u) |0) 

3 3 
+ [Xi{ui)^nj{u)aj4{u,ui) + X2{ui)Y^^lj{u)aj5{u,ui)]B3{u)\Qi) (3.32) 

j=2 j=2 
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So, will be an eigenstate of t{u) with eigenvalue 



3 

Ai(u, Ml) = ^'^flj{u)Xj{u)aji{u,ui) (3.33) 

provided the following equation is satisfied 

Xijui) _ Ifj=i^j{u)aj3{u,ui) _ Ej=2%WQj5Km) _ Q^^^^ ^^^^^ 

Explicit calculations of these expressions for the models considered in this paper will be presented in the 
next section. 

3.2 The two-peirticle state 

For m = 2 we have to seek eigenstates of t{u) in the form 

^2{ui,U2) = Bi{u,)B,{u2) \0)+B2{u,)T{u„U2) |0) (3.35) 

where r(ui,'U2) is an operator- valued function. Next we will use the condition that ^'2(^1,^2) must be 
normal-ordered to find r(ui,U2)- 

The first term of the right hand side of (3.35) has its normal-ordered form gives by the commutations 
relations: 

Bi{ui)Bi{u2) = w(mi,W2) [Bi{u2)Biiui) + Gd^{u2,Ui)B2{u2)Vi{ui) + Gd^{u2,Ui)B2{:LL2)'D2{ui)] 

-Gd^{ui,U2)B2{ui)Vi{u2) - Gd2{ui,U2)B2{ui)V2{u2) (3.36) 

where 

Xz{ui - U2)Xi{ui - U2) - Xq{ui - U2)y6{ui - U2) /o o'7^ 

Ul{Ui,U2} = r -, r {6.6() 

Xl{U\ - U2)X3[Ul - U2) 

Gdi(wi,W2) = — ^ 7 i (3.38) 

a;3(ui - U2) x\{2u2) 

n ( \ 3^6(^1+^2) „„s 
GdM\,U2) = 7 ■ r (3.39) 

X2\Ux + U2) 

Here we have used the following identities valid for both models, 

y6(-w) x-i{%L)xii{u) - x^(u)y^(%l) 



Xz{-u) X3{u)xi{u) - XG{u)ye{u) 

X2{2u) _ y5{u - v)x2{u + v) + X2{u - v)x5{u + v)fi{u) 
xi{2u) y5{u - v)xi{u + v) 



(3.40) 



(3.41) 
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Now we can see that (3.35) is normally ordered if it satisfies the condition 

*2(W2, Ml) = W(W2, Ui)*2(wi, U2) (3.42) 

This condition fixes r('Ui, U2) and, by construction, the unique candidate for eigenstate of t{u) in the case 
m = 2 has the form 

«'2(W1,U2) = Biiui)Biiu2) |0) + GdAui,U2)B2{ui)Vi{u2) \0) + Gd,{ui,U2)B2iui)V2iu2) |0) 

(3.43) 

The action of t{u) on this state in its normal ordered form is obtained with the following commutation 
relations, in addition to those presented in the m = 1 case: 

Viiu)B2iv) = bniu,v)B2iv)Vi{u) + bu{u,v)B2{u)Viiv) + b,:,iu,v)B2{u)V2{v) 

+bi4{u, v)B2{u)V3{v) + bi5{u, v)Bi{u)Bi{v) + bie{u, v)Bi{u)B3{v) (3.44) 

V2{u)B2iv) = b2i{u,v)B2{v)V2{u) + b22iu,v)B2iu)Vi{v) + b23iu,v)B2iu)V2{v) 

+b24{u,v)B2{u)V3{v) + b25{u,v)Bi{u)Bi{v) + b26{u,v)Bi{u)B3{v) 

+b27{u,v)B3{u)Bi{v) + b28{u,v)B3{u)B3iv) (3.45) 

V3{U)B2{V) = b3i{u,v)B2{v)V3{u) + b32{u,v)B2{u)Vi{v) + b33{u, v)B2{u)V2{v) 

+634 (w, v)B2{u)V3{v) + b35{u, v)Bi{v)Bi{v) + 636 (u, v)Bi{u)B3{v) 

+b37{u, v)B3iu)Bi{v) + bssiu, v)B3{u)B3{v) (3.46) 

Ct{u)Bt{v) = cu{u,v)Bi{v)Ct{u) + cr2{u,v)Br{v)C3{v) + ct3{u,v)Br{u)C3{v) 

+Ci4,{u,v)B2{u)C3{v) + C-[_r,[u,v)B2{v)C2{u) + ci&{u,v)Vi{v)Vi{u) 

+Ci7{u,v)T>i{v)T)2{u) + cis{u, v)T)i{u)T)i{v) + cig{u,v)Vi{u)V2{v) 
+ciio{u,v)V2{u)Vi{v) + cui{u,v)V2{u)V2{v) (3.47) 

C3{u)Bi{v) = C2i{u,v)Bi{v)Ci{u) + C22{u,v)Bi{v)C3{v) + C23{u,v)Bi{u)C3{v) 

+C24{u,v)B2{u)C3{v) + C25iu,v)B2iv)C2iu) + C26iu,v)'Di{v)Vi{u) 

+C27(m, v)T>i{v)V2iu) + C2s{u, (m)I?3 (w) + C29 (u, w)!?! (w)I?2 (^') 

+C2ioiu,v)Vi {u)V2 {v) + C2ll{u,v)V2{u)'Di{v) + C212 {u, v)T>2 {u)V2 [v) 
+C213(m, v)V3{u)V^{v) + C2ia{u, v)V3{u)Vi{v)) (3.48) 

After a straightforward calculation we obtain 

3 

t{u)'^2{ui,U2) = ^^^VLj{u)Xj{u)aji{u,Ui)aji{u,U2)]^2{ui,U2) 
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+[an{ui,U2)Xi{ui)'^Qj{u)aj2{u,ui) + a2i{ui,U2)X2{ui)'^Qj{u)aj3{u,ui)]Bi{u)Bi{^^^ |0) 
3=1 j=i 
3 3 

+[aii{ui,U2)Xi{ui)'^Qj{u)aj4{u,ui) + a2i{ui,U2)X2{ui)'^Qj{u)aj5{u,ui)]B3{u)Bi{^^^ |0) 

3=2 j=2 

3 3 
+ [aii{u2,Ui)Xi{u2)'^^j{u)aj2{u,U2) + a2l{U2,Ui)X2{u2)'^^lj{u)aj3{u,U2)]u!{ui,U2)Bi{u)Bi{ui) |0) 
3=1 3=1 
3 3 
+ [ail{u2, Ui)Xi{u2) ^ il.j{u)aj4{u, U2) + 021(^2, Ui)X2{u2) ^ f2j(M)aj-5(M, U2)\ljj{ui,U2)Bz{u)Bi{ui) |0) 
j=2 j=2 

3 3 
+ [Xx{u{)Xi{u2)^^j{u)Hji{ui,U2) + A'i(ui)Ai(u2)^Oj(u)iIj3(«i,M2) 
j=l 3=1 

3 3 

+A:'2(Ml)Afi(«2)Xl^j(")^i2(ui,U2) + A'2(Mi)Af2(«2)Xl^iW^i4(wi,U2)]fi2(^^ |0) (3.49) 
j=l 3=1 

where 

ifll(-«l,U2) = ai4(u, Ml) (ci6(mi,M2) + Ci8(mi,U2)) + ai5(u, Ml) (C26(W1,U2) + C29(U1,U2)) 

+6i2(m, ui)Gdi (wi, M2) + M)aii(u, Mi)ai2(u, ■U2)Gdi (u, w-i) 

i?12(Wl, W2) = ai4{u,Ui) (ci7(mi,M2) + Ciio(mi,M2)) + ai5(u,Mi) (c27(mi,U2) + C21i(mi,M2)) 

+6i3(u, wi)Gdi (ui, U2) + u)aii (u, wi)ai2(M, U2)Gd^ (u, ui) 

-H'lsWl, W2) = ai4(w, Ui)ci9(ui, U2) + dlsCw, Ml)C21o(wi, W2) + 612 (U, Ui)Gd2(wi) '"2) 

+uj{ui,u)a\i{u, ui)ais{u, U2)Gd^ {u, ui) 

Hi4.{ui,U2) = aii{u,U\)c\ii{u\,U2) + ai5(u,Ui)C212('Ul,W2) + h\z{u,Ui)Gd^{u\,U2) 

+Lo{uuu)aii{u, ui)ais{u, U2)Gd2{u, ui) (3.50) 
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and 

Hji{ui,U2) = aj6{u,Ui) (C16(U1,M2) + Cis{ui,U2)) + aj7{u,Ui) {C26{U1,U2) + 029(^1, M2)) 

+bj2{u, (wi, «2) + u)aji {u, Ui)aj2{u, U2)Gdi {u, ui) 

+aji{u, ui)aj4{u, U2)di3{ui,u) 

Hj2{ui,U2) = aj(i{u,ui) {cn{ui,U2) + ciw{ui,U2)) + ajr{u,ui) (027(^1, U2) + C2ii(mi, 'it2)) 
+bj3{u, uijGdi U2) + uj{ui,u)aji{u, Ui)aj2{u, U2)Gd2 {u, ui) 
+aji{u, ui)a,ji{u, U2)du{ui,u) 

Hj3{ui,U2) = aj6{u,Ui)cig{ui,U2) + OjT (u, Ml )C210 («! , M2) + 6j2 (u, Ml )G<i2 (ui , U2) 

+u){ui, u)aji{u, ui)aj3{u, U2)Gdi {u, ui) + aji{u, Ui)aj5{u, U2)di3{ui,u) 

Hji{u\,U2) = aj6{u,Ui)cin{ui,U2) + aj7{u,Ui)C212{ui,U2) + bj3{u,Ui)Gd2{ui,U2) 

+u){ui, u)aji {u, ui)aj3{u, U2)Gd2 («, wi) + aji(u, Mi)aj5(u, U2)du{ui,u) (3.51) 

for j = 2, 3. 

Again, '^2{ui,U2) will be an eigenstate oit{u) with eigenvalue 

3 

A2(u, ui,U2) = ^r2j(u)A'j(u)aji(u, ui)aji(u,U2) (3.52) 

provided the following equations are satisfied 

Xi{u{) ~.a2i{ui,U2) Xi{u2) sa2i{u2,u-i) 

' ' = e{u,ui) -, —j-^ = Q{u,u2) 3.53) 

X2{Ui) au{Ui,U2) X2{U2) aii{U2,Ui) 

where 0(ui), i = 1, 2 are given by (3.34). 
3.3 The n-particle state 

Prom the previous results we can are seek for operator valued functions with a recurrence relation of the 
form 

^n{u,...,Un) = Bi{ui)^n-l{U2,--- ,Un) 

n 

+S2(ui) ^ri'^(Wl, . . . ,M„)$„_2(W2, • • • ,Ui,... ,Un)Vi{ui) 



+S2(ui)^r2'''(Ml,-- - ,U«)*n-2(w2, • • • ,Un)'D2{Ui) (3.54) 

i=2 

It was shown in [6] that the above operator is normal ordered satisfying n — 1 exchange conditions 

$„(Ul, ... ,Ui, Ui+i, ... ,Un) = Oj{Ui, Mi+i)$„(wi, . . . , Wj+l, Uj, ... , M„) (3.55) 
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provided the functions Fa\ui, ... , u„) are given by 



i-i 



ri'-'(Mi,... ,u„) = Yl (^ci{ui,Uk)Gd^{ui,Ui), (a =1,2) (3.56) 

j=2 k=2,kjti 

Therefore the n-particle state will be given by 

*„(mi, ... ,Un) =<^„{u,... , u„) |0) (3.57) 
and the action of the operators Va{u), a = 1, 2, 3, on this state will be represented by 

n 

X>a(w)*n(wi, • • • ,Un) = A'a (u) ttal (u, Wi)*n(Ml, • • • ,Un) 

i=l 

n i—1 n n 

+ ^ JJ w(uj,«i)[A'i(u)a„2(w,Wi) JJaii(u,Uj) + A'2(u)aa3(w,Wj) JJa2i(u,Uj)]Bi(u)*„_i(uj) 

1=1 j=l jT^i jjii 

n i~l n n 

+ (1 - (5a, l) ^ Ul{Uj,Ui)[X-i{u)aa4{u, Ui) aii(u, Uj) + X2{u)aa5{u, Ui) JJ a2l(w, Uj)]B3(u)*„_i(Uj) 
i=l j=l jT^i jT^i 

n—1 n I n n 

+ < A'i(uj)A'i(uj) JJ an{ui,Uk) JJ an{uj,ui)Ho:i{ui,Uj) 

i=l j=i+l ^ fcT^iJ ^T^iJ 

n n 
+A'2(Wi)A'i(Uj) JJ a2l(Wi,'Ufe) au{Uj,Ul)Ha2{Ui,Uj) 

n n 
+Xi{Ui)X2{Uj) JJ ail(Ui,'Ufe) Y\_ 0,2l{Uj,Ul)Ha3{Ui,Uj) 
k^i,i IjtiJ 
n n 
+X2{Ui)X2{Uj) JJ a2l{Ui,Uk) J\ a2l{Uj,Ul)Ha4{Ui,Uj 
k^i,j 
i-1 i-1 

xJJa;(ufe,Ui) JJ u{ui,Uj)B2{u)^n-2{ui,Uj) (3.58) 

fe=l l=ljii 

Finally, the corresponding n-particle eigenvalue problem will be 

t{u)^n{ui,-- - ,Un) = i^^aiu)Xa{u)Y\_aaliu,Ui)\'^n{ui,... ,Un) (3.59) 
\a=l i=l / 

provided that the BA equations are satisfied 

— — - = e{uk) II r, (fc = l,2,... ,n) (3.60) 

-^2 K j_^j^i^aii{uk,Uj) 
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4 Explicit solutions 



In this section explicit expressions of the eigenvalue problem are presented for both models,. First we 
recall the fundamental relation (3.28) to get the coefficients aij{u,v) which appear effectively in the BA 
expressions (3.59) and (3.60): 

Xi{u — v) X2{u + v) 



aii{u,v) 



a2i{u,v) = -ui{u,v)[ 



X2{u — v) xi{u + v) 

xi{u + v)xi{u + v)+ x^{u + v)y^{u + v) 



Xi{u + v)X2{u + v) 



X2{u - v) X2{u + vY + Xfi{u + v)yQ{u + v) 
a?,i{u,v) = — — — — - — ■ — J (4.1) 

Xs{U — V) X2(U + V)X3{U + V) 

For the factor with the boundary contributions (3.34) wc will consider only the expression 

Q/ X ^ ^2(u)a25{u,Ui} + n3{u)a35{u,Ui) 
* i}2{u)a24{u,Ui) + n3{u)a34{u,Ui) 

where the ^j{u) are given by (3.21) and 

xe{u - v) X3{u + v) xe{u + v) 

x3(u — V) X2{u + V) X2{U + V) 

xq{u + v) 



a25{u,v) 



a34{u,v) = -f3{u)[fiiv) 



X2{U + V) 

Xe{u + v) X6{u - v) X3{u + v) 



X2{u + V) X3{U - v) X2{u + V) 

y6{u - v) .xe{u + v)ye{u + v)+ xl{u + v). 



X3{u-V) X2{u + v)X3{u + v) 

xeju - v)y(i{u - v) + xl{u - v) ye{u + v) 

xl{u — v) X2{u + v) 

X6{u + v) yeiu - v) X6{u + v)y6{u + v) + xl{u + v) 

X2{U + V) X3{U-V) X2{U + V)X3{U + V) 

with the fi{u) given by (3.13). 

4.1 sl(2|l)(2) model 

Substituting the matrix elements of the TZ matrix (2.5) and the matrix elements of the K matrices 
(2.7) and (2.9) for this model we get the following expressions in our algebraic BA: first the coefficients 
aji{u, Ui) for the eigenvalue A„(u, {uj}) 

sinh(u + ?ii) sinh(?i — Ua — 2r]) 

an {u, Ui) = . — — — ^ , (4.4 

smh(M + Ui + 2r]) smh(w — Ui) 

sinh(u + Ui) sinh(u — Ui — 2ry) cosh(u + Ui + 3r]) cosh(u — ut + rj) 

a2l{u, Ui) = . , , r-TT ^ TT ^ TT 7) (4-5) 

smh(u + Ui + 2r/) smh(u — Ui) cosh(u + Ui + r]) cosh(u — Ui —r]) 
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, , cosh(M + Ui + 3ri) cosh(u — m + n) 

a3i{u,Ui) = -y— -— ^ 

cosn(u + Ui+'r]) cosn(u — Ui — r]) 

Second, the factors of A„(u, {uj}) and of the BA equations with boundary contributions 

^ , , cosh(2M + Sr?) a sinh u — 2 cosh u a cosh(M, + 77) — 2 sinh('u + 77) 

cosh(2u + 77) a sinh(w + 2ri) — 2 cosh(w + 2ri) a cosh(w + ry) + 2 sinh(u + ry) ' 

^ siiih(2(/ + 2?7) a siiih (/ — 2 cosh ('4 7') 

^^"'^ ~ sinh(2u) a sinh(u + 2??) - 2 cosh(u + 2j?) ' ^ ' ' 

a sinh u — 2 cosh u 

^^"-^ " ~ a sinh(w + 2??) - 2 cosh(u + 277) ' ^ ' 



P sinh u + 2 cosh m ^ (u) 
/3 sinh u — 2 cosh u f^{u) 



sinh(2M) /3sinh(u + 277) - 2cosh(u + 2j7) a::|^(u) 
" sinh(2u + 277) /3sinhw-2coshu /^(u) ' ^ ' 



^ cosh(2(/ — ;/) ,isiuh(;; + 2//) — 2c:orih((/ + 2;/) ,icosh((/ + ;/) — 2riiuh((; + tj) (u) 

^ cosh(2u + 77) /3 sinh u — 2 cosh u /3 cosh(u — 77) + 2 sinh(u — 7/) f^{u)' 

and 

_ sinh(2Mi + 27;) a cosh(Mi + 77) + 2 sinh(M» + 77) (412) 
* sinh(2ui) a cosh(ui + 77) — 2 sinh(ui + 77) 

From these data one can see that the n-particle state ^n{{ui}) is an eigenfuction of the transfer 

matrix (3.20) for the sZ(2|l)(^) vertex model with eigenvalue 

-A- sinh(u + Mi - 77) sinh(u - - 77) 
fj^ smh(u + Mi + 77) smh(u — Uj + 77) 

o / w f \TT sinh(u + - 77) sinh(u - m - T]) 
sinn(M + ttj + 77) sinh(u — + 77) 

cosh(w, + Ui + 277) cosh(u — Ui + 2r]) 



X 



cosh(u + ?ii) cosh(u — Mi) 



+ns{u)Mu) IT """"^It ^f^^ '""^^17 (4.13) 

cosh(u + Ui) cosh(u — m) 



provided that the parameters {ui} satisfy the BA equations 



sinh(uj + 77 \ ^''^ a cosh Ui + 2 sinh Ui (3 sinh(Mj + 7;) — 2 cosh(u, + 77) 



sinh('Ui —T)) J cosh u, — 2 sinh m (3 sinh(uj — 77) — 2 cosh(uj — 77) 

n" COrih((/, + ILj + 77) cosh((/,; — Uj + 77) 
^ . J ^ cosh(ui + Uj — rj) cosh(ui — Uj — rf) 

i = 1,2,..., n (4.14) 
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where wc have used the shifts Ui ^ m = Ui ^ rj to bring these expressions into a symmetric form in rj. 

The formulation of this model in terms of the QISM presented here is new. However, one can verify 
that our results give the energy eigenspectrum previously obtained in the framework of coordinate BA 
by Fireman et al. [18]. 

4.2 osp(2|l) model 

For this model the K matrices have no free parameters but we have to consider three cases. For both 
cases the coefficients aji{u,Ui) are given by 

sin(M + Ui) sinh(u — Ui — 2ri) 



aii{u,Ui) 
a2i{u,Ui) 



sinh(u + Ui + 2rj) sinh(u — m) ' 

sin(i;, + Ui + 4?]) sinh(u + 1*^ + 77) sin(u — + 2rf) sinh(u — Ui — rf) 
sinh(M + Ui + 3ri) sinh(M + Ui + 21]) sinh(u — Ui) sinh(M — Uj + ry) ' 



, , sinh(w + Mi + 5?7) sinh(M - + 377) /.irN 

a3i{u,Ui) = -—- — ; riT^^-r? \ — r- (4-15) 

smh(u + Ui + 3r]) smh('U — ut + t]) 

Though, the boundary contributions are different for each cases: 
4.2.1 The (1,M) case 



In this case we have 



sinh(2u + r]) sinh(2M + 677) 
sinh(2w + 2r]) sinh(2u + 877) 



= _^., sinh(2. + 6,) 
^ ' smh(2u + 4?7) 

03(u) = e^" (4.16) 



_ 2,, sinh(2») xl^ju) 
' sinh(2u + 2r/) f«{v) 



sinh(2(/) siuh(2(; + 0//) .xg^ (u) 
sinh(2u + 47?) sinh(2u + 3r?) /^(u) 



and 



9(..)=e^'' ""'^^;,^f^\ ^ = l,...,n (4.18) 
smh(2ui) 

Therefore, the 7i-particle state ^n{{ui}) is an eigenfuction of the transfer matrix (3.20) for the 
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osp{2\l) vertex model with boundaries (1,M). The corresponding eigenvalue is given by 

sinh(2w + r]) sinh(2w, + Grj) (u) -A- sin(w, + Ui) sinh(u — Ui — 2r]) 
n{u,{ui}) = sij^jj(2?i + 2r]) sinh(2M + 877) f^{u) jLl sinh(M + w., + 2??) sinh(u - Ui) 

sinh(2M + 677) sinh(2M) X2^{u)-rj sin(M + + 477) sinh(u + + 77) 



n 



sinh(2u + 477) sinh(2M + 2?7) f^{u) sinh(u + Ui + 877) sinh(u + Mi + 2r;) 

sin(u — Ui + 2ri) sinh(u — Ui — tj) 
sinh('U — Ui) sinh('U — Uj + 77) 

sinh(2u) sinh(2M + 5?7) xl^{u) -A- sinh(7i + 77,,; + 577) sinh(u — Ui + 877) 



n 



sinh(2u + 477) sinh(2u + 877) f^{u) sinh(u + Ui + 877) sinh(u — Uj + 77) 

(4.19) 



provided that its parameters {uijare solutions of the BA equations 



sinh(Mi + 2r;) \ '^^ -pr sin(7i,; + uj + irj) sinh(?ii + Uj + 77) sin(7ii — uj + 2ri) sinh(7ii — Uj — 77) 

sinhuj / { '^\ ^ sinh(iti + Tij + 8r;) sin(iti + 7ij) sinh(iti — lij — 2r;) sinh(uj — Mj + 77) 

i = l,..,n (4.20) 
4.2.2 The (F+,G+) case 



Here we have 

rii(tt) 

02(u) 



2„sinh(2M + 677) sinh(u + |r7) cosh(u + 577) 
sinh(2u + 277) sinh(u + §77) cosh(u + §77) 

sinh(27x + 6r/) sinh(u + fr;) 
sinh(27x + 4r/) sinh(u + ^rj) 



^sinh(7i + |r/) 
sinh(u + §77) 



^,{u) = _e-2«-4.=l^ (4.21) 



- ,-2u Sinh(M+|77) xf'iu) 
^^^> ~ sinh(«-|77) /^(u) 

sinh(2M) sinh(u + \ri) xf^ {u) 



■^^^"^ sinh(2w + 277) sinh(w - §77) /^(u) 

X iu) = ^,2u+4r, sinh(2M) sinh(M + ^77) cosh(77 + §7;) xf^ [u] 

' sinh(27x + 477)sinh(M-|7;)cosh(w+|7/) /^(ti) ^' ' 



and 



efu-) = sinh(2ui + 277) sinh(M, + \-q) 

sinh(27Xi) sinh(uj + §77) 



(4.28) 
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Therefore, the n-particle state "^ndui}) is an cigenfuction of the transfer matrix (3.20) for the osp{2\l) 
vertex model with boundaries {F^ , G+). The corresponding eigenvalue is given by 

smh(?i + Ui + r]) smh(u — Ui + rj) 

,0 t w t ^ TT + w» + 3??) sin(M -Ui + Srj) 

smh(u + Ui + 2r]) smh(u — Ui + r]) 

sinh(u + Ui) sinh(u — Uj) 



sinh(u + Ui + r]) sinh('U — Uj + 2r]) 



+ns{u)X,{u) fr «inh(H + ..+4,)sinh(.-.. + 4r;) 

with the BA equations 

'sinh(uj +77)^''^ /sinh(uj - i?7) \ ^ -A- sin(ui + + 27?) sin(u, - ■ + 2?]) 



n 



sinh(ui - r?) / Vsinh(Mi -\-\r])) .4.f , sinh(ui + Uj + r\) sinh(ui - Uj - 2r]) 

sinh(w, + Uj — rj) sinh((/,; — Uj — 77) 
sin(ui + Uj + rj) sinh(ui — Uj + rj) 
i = l,...,n (4.25) 

Again, A„('U, {«,}) and the Bethe equations have been written in their symmetric form. 

4.2.3 The (F-,G-) case 



Here we have 



, sinh(2u + 677) cosh(u + ^rj) sinh(u + ^rj) 



= , ...... . 9 , . ,3 



sinh(2u + 277) cosh(u + |7?) sinh(u + |r/) 
sinh(2u + 677) cosh(?i + |?7) 



sinh(2?/ + 477) cosh(M + ^r]) 



_2„COSh(w+|7?) xl^{u) 



Xi{u) = e 



cosh(M — I??) f^{u) 
sinh(2M) cosh(7i + ^77) x'^'^ (u) 



sinh(2u + 277) cosh(M - §77) /^(m) 



' sinh(2u + 477) cosh(w - §77) sinh(M + §77) f^{u) ^ ' ' 

= ^-.„.sinh(2.. + 277)cosh(.. + |,) 
smh(2uj) cosh(uj + 5??) 
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Therefore, the n-particle state '^n{{ui}) is an cigenfuction of the transfer matrix (3.20) for the osp{2\l) 
vertex model with boundaries {F^ , G^)- The corresponding eigenvalue is given by 



An{u,{ui}) = Qi(u)A'i(u) Jl 



sin(u + Ui — rj) sinh(u — Ui — rj) 



■Mj- sinh(u + Ui + r]) sinh(u — Ui + rj) 

, r. r \-v I ^ TT sin(u + Mi + 3??) sin(u - + 37?) 

smh(M + Ui + zr]) smh(M — m + rj) 

sinh(u + Ui) sinh(u — Ui) 



sinh('U + Ui + ri) sinh('U — Uj + 2r]) 

I r. / ^ V / ^ TT ^'"^(" + + '^^) sinh(M - Ui + 4r]) 

+ih(u)A:3(u) I I . , , — ^ . , , -— (4.29) 

f-J^smh{u + Ui + 2r])smh{u-Ui + 2r]) ^ ' 

and the BA equations are now given by 

' smh{ui + ri)\^'^ f cosh{ui — ^r])\^ jj- sin(Mi + Uj + 2?7) sin{ui — Uj + 2ri) 

?/.-■ — -?/.-■ — ri\ 

(4.30) 



sinh(ui - ??)/ Vcosh(ui + |r?) / , jfi sinh(uj + Uj + ??) sinh(Mi - Uj - 2r]) 

sinh(u, + Uj — r}) sinh(«j — Uj — rj) 



sin(ui + Uj — 2rj) sinh('Uj — Uj + rj) 

Here, both A„('u, {u,}) and the BA equation were written with Ui ^ Ui — r]. 
These three cases were also considered in [18] via the coordinate BA. 

5 Conclusion 

Here, with the aid of previous works [6, 7], two of the three-state graded 19- vertex models have their 
boundary algebraic BA derived using a generalization of the Tarasov's approach [17]. From our results 
for the transfer matrix one can in principle derive the free-energy thermodynamics, the quasi-particle 
excitations behavior as well as the classes of universality governing the criticality of gapless regimes with 
integrable boundary conditions. Moreover, the rather universal formula we obtained for the eigenvectors 
could be useful in future computations of off-shell properties such as form factors and correlation functions 
with boundary conditions of relevant operators. 

The algebraic BA for n-state models with periodic boundary conditions was developed by Martins in 
[19]. In a recent paper Galleas and Martins [20] have presented the algebraic BA for the vertex models 
based on superalgebras. Therefore we believe that the Martins's approach can be generalized to include 
the diagonal open boundary conditions. 

Finally we observe that the vertex models discussed in this paper share a common algebraic struc- 
ture with the non-graded 19-vertex models, the Zamolodchikov-Fateev [21] and the Izergin-Korepin [22] 
models. Thus, we expected that a fusion procedure for the sZ(2|l)*^^) model as well as the analytical BA 
formulation based in the quantum group invariance of the osp{2\l) model can reproduce our results. 
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